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^ ■ Abstract 

By using the point canonical transformation approach in a manner distinct from 
previous ones, we generate some new exactly solvable or quasi-exactly solvable poten- 
tials for the one-dimensional Schrodinger equation with a position-dependent effective 
Qh! mass. In the latter case, SUSYQM techniques provide us with some additional new 
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1 Introduction 



In recent years, quantum mechanical systems with a position-dependent effective mass 
(PDEM) have attracted a lot of attention due to their relevance in describing the physics 
of many microstructures of current interest, such as compositionally graded crystals [1], 
quantum dots [2], 3 He clusters [3], quantum liquids [4], metal clusters [5], etc. 

As in the constant-mass case, exact solutions play an important role because they 
may provide both a conceptual understanding of some physical phenomena and a test- 
ing ground for some approximation schemes. Many recent developments have been devoted 
to constructing exactly solvable (ES), quasi-exactly solvable (QES) or conditionally-exactly 
solvable potentials for the PDEM Schrodinger equation [6]— [17] by using point canonical 
transformations (PCT), Lie algebraic techniques or supersymmetric quantum mechanical 
(SUSYQM) methods. 

In this Letter, we will show that new ES or QES potentials in a PDEM background 
may be generated by using the PCT approach in a manner distinct from previous ones. We 
will then combine such results with SUSYQM methods to produce some additional QES 
potentials. 

2 PCT approach in a PDEM context 

As is well known (see, e.g., [14]), the general Hermitian PDEM Hamiltonian, initially 
proposed by von Roos [18] in terms of three ambiguity parameters a, f3, 7 such that 
a + [3 + 7 = —1, gives rise to the (time-independent) Schrodinger equation 

d 1 d 



Hij)(x) = 
where the effective potential 



dx M(x) dx 



+ V eS (x) 



%jj{x) = Eijj(x), 



V cS (x) = V(x) + \{(3 + 1)^ - [a(a + (3 + 1) + (3 + 1]^ (2) 

depends on some mass terms. Here a prime denotes derivative with respect to x, M(x) is 
the dimensionless form of the mass function m(x) = m M(x) and we have set ft = 2m = 1. 



2 



Let us look for solutions of Eq. (1) of the form 

*P(x) = f(x)F(g(x)), (3) 

where f(x), g(x) are two so far undetermined functions and F(g) satisfies a second-order 
differential equation 

F + Q(g)F + R(g)F = 0, (4) 

where a dot denotes derivative with respect to g. Since in this Letter we shall be interested 
in bound-state wavef unctions, we shall actually restrict ourselves to polynomial solutions 
of Eq. (4). 

On inserting Eq. (3) in Eq. (1) and comparing the result with Eq. (4), we arrive at two 

expressions for Q(g(x)) and R(g(x)) in terms of E — V e g(x) and of M{x), f(x), g(x) and 

their derivatives. The former allows us to calculate f(x), which is given by 

[M\ 1/2 (l [9W \ 
f(x) oc I — I exp I - J Q(u) du I , (5) 

while the latter leads to the equation 

E - V ^ x) = 2M g -'-AM\- g 7) + m{ R -2 Q -A Q )'2M^ + UP- (6) 
It is clear that we need to find some functions M(x), g(x) ensuring the presence of a 
constant term on the right-hand side of Eq. (6) to compensate E on its left-hand side and 
giving rise to an effective potential V e g(x) with well-behaved wavef unctions. 

In the constant-mass case, i.e., for M(x) = 1, this procedure has been thoroughly 
investigated [19, 20]. A similar study in the PDEM context looks more involved for two 
reasons: (i) there are now two unknown functions instead of only one and (ii) the usual 
square-integrability condition for bound-state wavefunctions has to be completed by the 
additional restriction \ip(x)\ 2 / \J M(x) — > at the end points of the definition interval of 
V(x) to ensure the Hermiticity of H in the Hilbert space spanned by its eigenf unctions [17]. 

In most applications of PCT that have been carried out so far in the PDEM context, 



the choice M = \g' 2 or g(x) = (1/A) J x \JM(u) du + v (where A, v are some constants) 
has been made (see, e.g., [12, 13, 14]). In the next two sections, we will explore the new 
possibilities offered by two other options, namely M = Xg' and M = X/g' or, equivalently, 
g {x) = (1/A) J x M{u) du + u and g{x) = (1/A) j x [M(u)}' 1 du + v. 



3 Generation of ES potentials in the M = Xg' case 



Substituting M = Xg' into Eq. (6) leads to 



E-V eS (x) = ^g' (r-±Q-±Q* 



(7) 



Some simple and interesting results can be derived from this relation by assuming that F(g) 
is either a Jacobi or a generalized Laguerre polynomial [21]. 
For F n (g) oc P^' b \g), n = 0, 1, 2, . . . , a, b > -1, we obtain 



r-\q-\q 2 



n(n + a + b+l) 



+ 



1-g* 
9 1 



+ 



1 



1-9 2 ) 2 
(b-a)(b + a) - 



\{a + b + 2)-\{b-af 



9' 



1 



a + b)(a + b + 2). (8) 



:i-^)22 v /v ' (l-g 2 ) 2 ^ 

A constant term can therefore be generated on the right-hand side of Eq. (7) by assuming 
g'/[X(l — g 2 )} = C, where C must be restricted to positive values in order to get increasing 
energy eigenvalues for successive n values. The solution of this first-order differential equa- 
tion for g(x) leading to a positive mass function reads g(x) = tanhgx, where q = XC > 0. 
Without loss of generality, we may set C = q 2 so that X — l/q. Hence we get 



g(x) = tanhgx, M(x) = sech 2 qx, 



Equations (5), (7), (8) and (9) then yield 



-co < x < +oo. 



(9) 



E„ 



V cS (x) 



q [n + 



b\ ( a + b + 2\ T7 



^a 2 + b 2 ) 



cosh 2 qx + - (a — b) (a + b) sinh qx cosh grr j 



Vn 



[a 2 - l)e 2qx + (b 2 - l)e- 2qx + a 2 + b 2 - 2] + V , 
ip n (x) oc (1 -tanhgx) (a+1)/2 (l +twhqx) {b+1)/2 P^ b \tanhqx), 



(10) 

(11) 
(12) 



where n — 0, 1, 2, . . . , V denotes some constant and we have to assume a,b > —1/2 in 
order to satisfy the conditions on bound-state wavefunctions in a PDEM context (observe 
that the square-integrability condition alone does not impose any restriction on a, b\). 
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By proceeding similarly for F n (g) oc L$(g), n = 0, 1, 2, . . . , a > — 1, from the 
relation [21] 



1 ■ 1 2 2n + q + l ( a + l)(q-l) 1 

R ~2 Q ~l Q = 2g 4? 4 (L,) 

and the condition g' / (Xg) = C > 0, we obtain the results 



g {x) = e~ qx , M(x) = e~ qx , -oo < x < +oo, (14) 

where we have set C = q 2 (hence A = — 1/q) and where without loss of generality we may 
assume q > 0. Furthermore 

E n = q 2 (n+^±^+V , (15) 

V cS (x) = ^q 2 [(a 2 - \Y X + e- qx \ + V , (16) 

Mx) oc exp|-i[(a + l)gx + e -^]}4 a )(e^), (17) 

where the PDEM background imposes an additional restriction a > —1/2 on the wavefunc- 
tions again. 

Turning now to the initial potential V(x), we find from Eq. (2) that V(x) = V c g(x) + 
q 2 [f(a, (3) cosh 2 qx - g(a,{3)} and V(x) = V cS (x) + \q 2 f(a,p)e c ' x , with f(a,(3) = (2a + 
l)(2a + 2 / 9 + 2) — 2a, g(a, j3) = (2a + l) 2 +(3(4a + l), for the Jacobi and generalized Laguerre 
polynomials, respectively. Hence, in both cases, for the choice of ambiguity parameters 
made by BenDaniel and Duke (a = 0, (3 = —1) [22], there is no distinction between 
V(x) and V e s(x). Furthermore, when the Jacobi polynomials reduce to Legendre ones, 
i.e., for a = b = 0, and the ambiguity parameters are those selected by Zhu and Kroemer 
(a = —1/2, (3 = 0) [23], V(x) becomes a constant potential V . Our results (10) and 
(12) then describe the generation of an infinite number of bound states for a free-particle 
potential in a sech 2 -mass environment [14]. For nonvanishing a, b values, Eqs. (10)-(12) 
may therefore be seen as a generalization of this interesting property. 

4 Generation of QES potentials in the M = \/g' case 

Whenever M = X/g', Eq. (6) becomes 

E - V ^ = T-h + i x{ R --2 Q -4 Q2 )- < 18 > 
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In such a case, we shall take for F(g) some polynomials of nonhypergeometric type satisfying 
the equation 

F + ^^F+^F = 0, (19) 

gi gS 

where we assume a, b, c, £ real, a ^ 0, b ^ and £ > 0. As shown elsewhere [24], this second- 
order differential equation has /cth-degree polynomial solutions provided b = —k(a + k — 1) 
and there exist k + 1 such solutions F n (g), n = 0, 1, . . . , /c, associated with /c + 1 distinct 
values c„ of c, if a is appropriately chosen. 
Substituting 

1 ■ 1 2 (2A; + a-2)(2A; + a) , c n , a(a - 3)£ 2 a 2 £ 4 



4- 4^2 % 3 + 2# 4 V ^ 

in Eq. (18), we find a constant term on the right-hand side of the transformed equation by 
choosing g' 3 / (\g 3 ) = C. Then with C = q 2 and A = q > 0, we obtain 



g(x) = e qx , M(x) = e~ qx , -oo < x < +oo. 



(21) 



Hence 



E n = q 2 c n + V , 
V cS (x) = q 2 

ip n (x) oc exp 



\{2k + a - 2){2k + a)e qx - -a{a - 3)£ 2 e - qx + \ 2 ^e- 3qx 
4 2 4 



\{a-2)qx + \aee- 2qx 



F n (e qx ) , 



(22) 
+ V , (23) 

(24) 



where n — 0, 1, . . . , k. 

The functions (24) turn out to be physically acceptable as bound-state wavefunctions 
provided a is restricted to the range a < —2k + |. We conclude that for such values and for 
the PDEM given in (21), the effective potentials (23) corresponding to k — 1, 2, 3, ... , are 
QES with k + 1 known eigenvalues and eigenfunctions. For k = 1 and k = 2, for instance, 
we find 



Co = ±af, Fo(g) oc g>±£, if a < --, 



(25) 



and 



Co = ^A£, ci = 0, 



a 



oc g 2 



a + 1 



£ 2 , ^ = \j2a{2a + 3), 
6 



if a < 



a + 2 
5 



(26) 



respectively. Observe on these two examples that for the values taken by a, £ and g(x), 
ijj n (x) has n zeros on the real line, so that ipo(x) is the ground-state wavefunction, while 
ip n (x), n = 1, 2, . . . , k, correspond to the nth excited states. 

The results presented here could be easily extended to more general polynomials. For 
instance, if instead of g 2 — £ 2 in (19), we had considered (g — g 3 ){g — g^) with g 3 and g 4 real 
but #4 7^ — #3, we would have obtained effective potentials containing an additional term 
proportional to e~ 2qx . 

Finally, it should be noticed that the PDEM being the same as that chosen for general- 
ized Laguerre polynomials in Sec. 3, the relation between V(x) and V c s(x) is also similar. 

5 SUSYQM approach 

Let us consider the intertwining relationship r\H = Hii], where H is the Hamiltonian defined 
in Eq. (1), Hi has the same kinetic energy term but an associated effective potential Vi ie g(x), 
and i] is a first-order intertwining operator rj = A{x)-^- + B{x). As shown in [14], such a 
relationship leads to the restrictions A(x) = M -1 / 2 and 

+ V lMx) = V* + ™ + ^-%, (27) 

with e denoting some arbitrary constant. 

A solution for B(x), which at the same time ensures that rj annihilates the ground- 
state wavefunction of H, is provided by B{x) = —ip' / '(y/Mtp ) together with e = E . 
In this (PDEM-extended) unbroken SUSYQM framework [25], the eigenvalues of Hi are 
Ei >n = E n+ i, n = 0, 1, 2, . . . , with the corresponding wavefunctions given by ipi >n oc r/Vwi- 

For the wavefunctions considered in Eq. (3), ip' /ipo in general contains two terms: 
V'o/V'o — /'// + Fq. In the ES potential case reviewed in Sec. 3, however, the second 
term vanishes since F (g) = 1, so that we obtain simple results for B(x), namely 

B{x) = -q[(a — b) cosh qx + (a + b + 2) sinh qx] (28) 
2 

and 

B(x) = -q[(a + l)e qx/2 - e"^ 2 ] (29) 
7 



for the Jacobi and generalized Laguerre polynomials, respectively. Substituting such func- 
tions in (27), we arrive at SUSY partners Vi je g(x), which have the same shape as V e g(x) 
and differ only in the parameters (ai = a + 1, b\ = b + 1 , V 0; i = V and a x = a + 1, 
— V + \q 2 , respectively). We conclude that the potentials V e g(x) are shape invariant. 
The QES potential case reviewed in Sec. 4 looks more interesting because F (g) being 
now a /cth-degree polynomial in g, the second term in ip' /ipo does not vanish anymore. As 
a consequence, the functions B(x) and V^ e g{x) become /c-dependent and given by 



B(x) = q 
V lte s(x) = V eS - q 2 



-(a- 2)e qx ' 2 + -a^V 3 ^ 2 - e 3 ^ 2 ^ 
2 V ' 2 s F 



L l a _ £ e ^ + V 2 ^* + 3e 2 ^ + 2e 3 ^ ^ - ^ 



(30) 



•(31) 



2 V 2/ 2 F yF F 

The SUSY partners Vi >e g(x) therefore contain some terms which are rational functions in 
e qx . For k — 1 and /c = 2, for instance, we obtain 



V 1:CS (x) = q 2 



1, 



+ 



l)(a + l)e 9a; 
3£ 2 2£ 3 



e gx _|_ £ ( e <pr _|_ £)2 



-a 2 ee~ qx + -aYe- 3qx 
2 s 4 s 



(32) 



and 



V^i,eff(a:) = q 2 



-(a + l)(a + 3)e qx - -a 2 i 2 e~ qx + -a 2 ^e~ 3qx 
4 2 4 



+ 



< 2 



2 M 



4a^ 



(a + 2) 3 ~ ±v ^ (a + 2) 4 
6(a + 2)(a+ l)e ga; - (a + 6)A£ 



+V + q 



2 A£ 
a + 2' 



2qx 



a+2 

(3a + 4)e 9a: - A£ 



2 ' 



(33) 



respectively. Such effective potentials provide us with some new examples of QES potentials 
in a PDEM environment with k known eigenvalues and eigenfunctions. 



6 Conclusion 



In this Letter, we have investigated the problem of the one-dimensional Schrodinger equa- 
tion in a PDEM background from several viewpoints. By using first the PCT approach and 
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assuming a relation between the new variable g = g(x) and the mass M(x) that differs from 
the usual one, we have constructed some new ES or QES potentials. The former are asso- 
ciated with either Jacobi or generalized Laguerre polynomials, while the latter correspond 
to some /cth-degree polynomials of nonhypergeo metric type. 

We have then considered an equivalent intertwining-operator approach and shown that 
while our ES potentials are shape invariant, the SUSY partners of our QES potentials are 
new. In the latter case, iterating the procedure would lead us to a hierarchy of SUSY 
partners with an increasingly complicated form. 

The method described here could be used to generate other classes of masses and po- 
tentials providing exact solutions of the PDEM Schrodinger equation. 
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